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Abstract 

By using the localized character of canonical coherent states, we give a straightforward derivation 
of the Bargmann integral representation of Wigner function {W). A non- integral representation is 
presented in terms of a quadratic form W oc V^FV, where F is a self-adjoint matrix whose entries 
are tabulated functions and V is a vector depending in a simple recursive way on the derivatives 
of the Bargmann function. Such a representation may be of use in numerical computations. We 
discuss a relation involving the geometry of Wigner function and the spacial uncertainty of the 
coherent state basis we use to represent it. 

PACS numbers: 03.65.-w 
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I. INTRODUCTION 



The Heisenberg uncertainty principle assures that no function of the canonically conju- 
gated pair {q^p) can be defined in a way to be interpreted as a genuine probabihty density 
in phase-space. Despite this fact, many phase-space representations of quantum mechan- 
ics have been developed and demonstrated to be powerfuU tools in different fields, such as 
semi classical limit 1^3, 1^, li] and quantum optics j^, Q] . 

In 1932 Wigner 7|] introduced his famous quasi-probability function W, so named for 
its possible negativity in some regions of phase-space. For a pure ensemble with density 
operator p = \ip){'ip\, W is usually presented as an integral in configuration space: 

W^(g,p) = ^y" dy{q + y/2\^){^\q-y/2)e-'^^/\ (1) 

Among other appealing properties, the above function yields the correct marginal probabili- 
ties: integration of W in the variable p gives the position probability density, and conversely, 
integration in the g-axis leads to the probability density in momentum space. The Wigner 
function contains all the information on the ensemble, providing an alternative formalism for 
quantum mechanics. In the case of mixed ensembles it enables the calculation of quantum 
averages in a classical way. 

Considering that the simplest way to address quantum mechanics in phase-space is, ar- 
guably, by employing coherent states, it is natural to ask how the Wigner function is related 
to this representation. Not surprisingly, this connection has been partially studied by Cahill 
and Glauber Q], and Wiinsche jol. 1lo|. on very formal grounds. We shall be concerned 
with the slightly different, but equivalent question of how the Wigner function is related 
to the analytical representation associated to the coherent states (the so-called Bargmann 
representation [HI]). 

Our purpose in this work is, firstly, to use the fact that coherent states are maximally 
localized structures in phase-space to give a simple derivation of the Bargmann integral 



representation of Wigner function (section II 
this result is simpler than the previous ones 



I. We believe the method we use to obtain 



10|. It provides more physical insight 



and leads to the non-integral representation given in section IV in a natural way. These 

discrete representations have been often used in numerical evaluations of Wigner function in 
□ 

theoretical and experimental situations, e. g., in the reconstruction of motional states of 
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trapped atoms and ions Q, [l^. In section V we discuss a relation involving the geometry of 
Wigner function and the spacial uncertainty of the coherent state basis we use to represent 
it. Our final comments are outlined in section VI. 



II. PRELIMINARY DEFINITIONS 

Canonical coherent states [l^, [l^ are defined in terms of number states associated to 
quantum harmonic oscillators: 

|.') = e-'*^'/2^^|^)^ (2) 
Vn! 

where \/2 z' = {q' /b + ibp' /H) is a complex label, q' and p' are the mean values of the related 
quantum operators in the state [zQ,^ and 6 is a positive arbitrary constant. Such states 
constitute an over-complete basis [17,] of the Hilbert space, and therefore, can be used to 
express the resolution of unit: 

where d^z'/vr = dg'dp'/27r^. We recall that in the position representation we have 



(3) 



from which it is clear that the parameter b is related to the uncertainty in position of the state 
\z') {b = ^/2Aq'). This constant can be freely chosen. Note, however, that the minimum 
uncertainty relation Aq'Ap' = H/2 must be satisfied (implying b = ^/2H/ Ap'). 

III. INTEGRAL REPRESENTATION REVISITED 

In order to define the Wigner function as a phase-space integral we start by conveniently 
inserting coherent-state unit operators in expression ([1]) 

Wiq,p) = -^ [ ^^{z"\i;){i;\z') [ dy {q + y/2\z"){z'\q - y/2) e-'^y/\ (4) 
Znn J 71 n J 

Since {q + y/2\z"){z'\q — y/2) is a Gaussian function of y (see equation (I3l)) the integration 
in this variable readily gives 

^-2z*z f j2 / f^]2 /> . . ^ 

W{q,p) = W{Z, Z*) = —— \ ^ ^ ^^//|^^^^|^/^g2.'*.+2.-."-.-."/2-.-.'/2-.'*." _ (5) 

Tifi Inn 



In order to better explore the properties of coherent-states we change variables as follows: 
z" = w and z' = w + 6w. Since the Jacobian determinant is unitary one gets 

W{z,z*) = — / ^ — - {w\i!){ilj\w + 6w) 

nh J n TT 

exp{-2w*w + 2z*w + 2zw* - Sw*5w/2 - {3w/2 - 2z)5w* - w*Sw/2} . (6) 

Note that, due to the minimum uncertainty character of coherent states, the integration in 
w and w* involving {w\ip) {iplw + 6w) , with ip being a function, falls-off in a Gaussian-like 
way for 6w 7^ 0. Therefore, it is safe to make an expansion around this point. However, 
before proceding to such an expansion, we recall that {iplw + 6w) is not an analytical 
function of w + 6w, since there is a (trivial) dependence on w* + 6w*, so that {ip\w + 6w) = 
fr(,{w + 6w,w* + 6w*). Thus, the referred expansion, in the coherent state representation, 
amounts for a two variable Taylor series. This unnecessary complication can be avoided 
in the Bargmann representation, defined as f{w) = exp{w*w/2){tlj\w), which provides a 



description based on entire functions Re- writing the Wigner function in this formalism 

we get 

-2z*2 r j2 A'^f^yA 
' Tlh J 71 TT 

(7) 

Since f{w + 5w) is an analytic function of its argument one can proceed a single variable 
Taylor expansion: f{w + 6w) = J27=o ^/(^) We obtain 

— 2z*Z p J2 ^ 1 

W{z,z*) = \ ^r{w)e-'-'-^'^'-^'^-*Y.-^^nf(^) 

n=0 

d [bw) ^^n^—Sw*5w—2{w—z)5w*~w*Sw 



TT 



The integral in Sw and Sw* is given by 

dw* 



thus 

wiz, = ^ / ^ r (-) e-*--^^^*- -)"d^/(-) ■ (10) 

By writing f{w) = YlT=o^kw'' we have ^/(w) = Efclo '^fc (fSt)!^^"''' where we used the 
fact that l/|r(— A^)| ^ for = 0, 1, 2, ... to extend the second sum in k from (n, oo) to 
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(0, oo). This leads to 




n\{k — n)\ 



k\ 



{2z - 2wYw 



OO 



k=0 



where we used the binomial expansion. We finally get 



W{z,z*) 




(11) 



which can be written in a more symmetrical way as: 



W{z, z*) 




nz + w/2)f{z-w/2) e 




(12) 



where f*{z + w/2) = exp{{z + w/2){z* + w* /2) /2}{z + w /2\ij) and f{z-w/2) = exp{{z - 
w/2){z* — w* /2)/2}{ip\z — w/2). This expression is the phase-space analogous of equation 
([T]). It is clear that in the case of a mixed ensemble of states {ipi) with statistical weights pi 
(p = '^iPi\4'i){4'i\) we have W = "^iPiWi, where each Wi is given by equation (fT2il 

IV. NON-INTEGRAL REPRESENTATION 

It is well-known that integrals over phase-space are, in some situations, not suitable 
for numerical evaluations of Wigner functions [12|]. In this section we obtain a non-integral 
representation of W, in terms of derivatives of Bargmann functions, from expansions around 
w = 0. The same argumentation used in the previous section is valid here: f*{z + w/2)f{z — 
w/2) falls-off in a Gaussian-like way for w ^ 0. Let us write 




(13) 



Substituting in (JT2l) we have 




(14) 



where 




(15) 
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The integration gives 



Qn 



4(_l)i2"+%!j!2*Ve- 



a=0 



dz 

-ir 



l-2a 



a\{n — cr)!(j — a)\ 



= 4(-l)^2"+^n!j!z*Ve-^*^ 2Foi-n, -j, -\z\-^) , 
where, again, the summation in a may be extended to oo by means of the already referred 



property of the Gamma function and, 2-^0 denotes the confluent Hypergeometric function of 

;t expressio] 



second kind 3, Going back to ffT^ we obtain the compact expression 



-2z*z 



W{z,z* 



J2 2Fo(-n,-j,-|zr^)^*V 



dz*"- dz^ 



(16) 



As a simple application we take a harmonic oscillator eigenstate {ip) ■ 
f{z) = z^/^fW\, leading to 

-2-*^ .M-n,-3,-\z\-^) 



\N). We then have 



W{z,z*) = N\ \z 



2N 



(N-j)\{N-n)< 



Summation in n and j readily gives the expected result 



W{z,z*) = N\ \2z 



2N 



-22*2 



TT 



{-iy{4:z*zy 

a\{N - a)\{N - a)\ 



\N 



-22*2 



-Ln{^z*z) 



where L^r denotes the Laguerre polynomial of degree N . 

We note that equation ffTUl) can be written more elegantly as the quadratic form 

„-22*2 



W{z,z*) 



Tltl 



V^FV 



(17) 



with the vector V and the Hermitian matrix F given by 



V 





/ 


\ 






d/ 








dz 
















dz^ 






v 




) 





2Fo(0,0) 2i^o(0,-l)z 2i^o(0,-2)^2 

2Fo(-l,0)z* 2F^{-l,-l)z*z 2Fo(-l,-2)^V 
2Fo(-2,0)z*2 2Fo(-2,-l) A 2i^o(-2,-2)z*V 



V 



(18) 



where we have suppressed the third argument (— |z|~^) in the matrix elements. Although 
expression (1171) does not seem to be particularly usefull in simple analytical calculations, 
it has potential value from both, formal and computational points of view. Equation ffT7|) 



may be suitable for numerical evaluations of W in situations where the Bargmann function 
and its derivatives are easier to obtain than the equivalent integrations in phase space or 
configuration space (which seems to be usually true). The matrix F is built with tabulated 
functions and the vector V has a simple recursive nature, namely Vj+i = dVj/dz. The 
diagonal elements of F can be put in a simpler form: F„_„ = Ln{z*z), which makes 
clear that the previous definitions of F and V are less effective in regions of large \z\. The 
alternative definitions = (/, zdf/dz, z^ d? f /dz"^ , ...) and F„j = 2FQ{—n,~j,—\z\~'^) 
(orthogonal matrix) with W = ^-^^V^FV may be more convenient, depending on the 
system under study, for regions of large \z\, where the entries of F and the derivatives of the 
Bargmann function are both small. 



V. W DOES NOT DEPEND ON b: GEOMETRICAL MEANING 

A final note on the use of variables z and z* to express W is in order. Since in the original 
definition ([T]) there is no mention to the parameter b, W must not depend on it regardless 
the variables we are using to cover the phase-space. Thus, we have dW/db = 0, which leads 

to 

db dz dz* ' 

where we have used dz/db = —z*/b and dz*/db = —z/b. It is clear that some explicit 
dependence on b is required in order to compensate the implicit dependence contained in 
z and z*, so that dW/db ^ 0. As we shall see, the above relation contains a piece of 
geometrical information. By writing 



— - — (b— - i-—^ —- — (b— t-—^ 
dz a/2 \ dq b dp J ^ dz* ^/2 \ dq b dp J ^ 

it is easy to show that both sides of equation f|T9l) are independent of b. We obtain 



(20) 



^dW _ dW _ dW 
db ^ dq ^ dp ' 



This relation can be put in a more suggestive form 

b dW 



n . grad W , (22) 



v/g2 + p2 db 

where n is the unit vector in phase-space parallel to {q,—p) and grad = {d/dq, d/dp) is 
the gradient operator in phase-space. We see, therefore, that bdW/db is proportional to 



the component of the gradient of the Wigner function in the direction (g, — p), i. e., to the 
variation of W in the n direction (see figure 1). 

As an illustration let us assume that is itself a coherent state \U), with \/2U = 
{Q/B + iBP/H). By using ( fT2l) (note that b and B not necessarily coincide) it is easy to 
obtain 

1 (B^-b\ 2 B^ + b" 



Wiz, z*) = - exp I ^^(^ +zn- 

B^ — B^ + b^ 

+ {zU + z*U*) + {zU* + z*U) - 2U*U 

Bb Bb 

It is clear from the corresponding expression in the {q,p) variables, W{q,p) = ^exp{ — (g- 

QY/B'^ — B'^ij) — P)'^/fi^}, that W is independent of b. In this particular case, relation ([22 

estates that bdW/db = 2W \-^{q - Q) + ^{p - P)] = (g, -p) . grad W , for any b. 



VI. CONCLUDING REMARKS 



In this work we presented a simple and intuitive derivation of the Bargmann integral 
representation of Wigner function along with a non-integral representation. We hope that 
expression (ITTl) will be of use in numerical computations for a class of theoretical and exper- 
imental problems in which the derivatives of the Bargmann functions are readily accessible 
or, at least, easier to calculate than the integrations in configuration space, equation ([T]), or 
phase-space, equation ffT^ . Finally we call attention to the fact that, despite the fact that 
W does not have an overall dependence on the parameter b, its explicit dependence has a 
simple geometrical meaning. 
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